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ABSTRACT 

We use archival data of NASA's Chandra X-ray telescope to compile an X-ray light 
curve of all four images of the quadruply lensed quasar Q2237-I-0305 (2q=1.695) from 
January 2006 to January 2007. We fit simulated point spread functions to the four 
individual quasar images using Cash's C-statistic to account for the Poisson nature 
of the X-ray signal. The quasar images display strong flux variations up to a factor 
of ~ 4 within one month. We can disentangle the intrinsic quasar variability from 
flux variations due to gravitational microlensing by looking at the flux ratios of the 
individual quasar images. Doing this, we find evidence for microlensing in image A. 
In particular, the time-sequence of the flux ratio A/B in the X-ray regime correlates 
with the corresponding sequence in the optical monitoring by OGLE in the V-band. 
The amplitudes in the X-ray light curve are larger. For the most prominent peak, the 
increase of the X-ray ratio A/B is larger by a factor ~ 1.6 compared to the signal 
in the optical. In agreement with theory and other observations of multiply imaged 
quasars, this suggests that the X-ray emission region of this quasar is significantly 
smaller than the optical emission region. 

Key words: gravitational lensing: micro - X-rays: galaxies - galaxies: quasars: indi- 
vidual Q2237-I-0305 - cosmology: observation. 



1 INTRODUCTION 



The 



. Q223 7-I-0305 was discovered in 1984 

(|Huchra et al.l 1 19851 ) during a spectroscopic survey of 
nearby galaxies. The spectrum of the nucleus of the barred 
spiral galaxy 2237+0305 (2g=0.0394) was found to be 
superimposed by a quasar component at a redshift of zq = 
1.695. 

The first high resolution obser vations of the system showe d 
three images of the quasar (|Tvson fc GorensteinI Il985l ). 
This number was soon corrected by l|l988l ) who first 
observed all four known point like quasar images around 
the core of the galaxy. The images are arranged in a nearly 
symmetric way, hence the name 'The Einstein Cross'. 
Spectroscopy proved t h at they are images of a single quasar 
jSchneider et al.lll988l : lAdam et al.lll989l ). The images are 
separated by up to 1.8" from each other, and are labeled 
A through D. 

After the prediction by IChang fc Refsdal J 19791 ) and 
initial theoretical studi es (e.g., Kayscr et_ 



microlensing' was the first evidence for microlensing 
generally, including stellar microlensing in the galaxy 
which was not discov e red until 1993 (lAubourg et al.l 
lAlcock et ai] Il993l : lUdalski et all Il993l ). Today 'The 
Einstein Cross' is one of the best studied multiply lensed 
quasars and ther e have been a lot of monitoring programs 
in the optical (ICo rrigan et al.' '1991'; 'Pen fc et al l Tl993l: 
0stensen et al. 1996: Wozniak ct al. 2000b a; Schmid t et al.l 



20021 : lAlcalde et al.lf2002l : lUdalski et al.l 1200 6) where much 
microlensing activity has be en observed. The q uasar has 

'l998f). the NIR 



also been detecte d in the U V jBlanton et aL 
lAgol et al.ll2009l). the MI R (|Agol et al.ll200oi ) and the radio 



i (1 19791) 
ct afi |j 



. . . _ 19861 ): 

[Schneider fc Weisd (|l986l ): IPaczynskil (1 19861)). microlens- 
ing was first detected i n Q2237-I-0305 (|lrwin et"all 1 19891 : 
IWambsganss et al.|[l990l ). In fact, this detection of 'quasar 



regime ( Falco et al.lll99^ ). However, there is no published 
light curve of Q2237+0305 in the X-ray regime yet. 
The X-ray emission of Q2237-I-0305 was first detected with 
ROS AT/HRI observations in 1997 (|Wambsganss et all 
Il999h . Since then there have been other X-ray observations 
of Q2237+0305 like a single spectroscopic observation 
with XMM-Ne wton from 2002 (data set ID: 0110960101; 
PI: Watson) (iFedorova et al.l |2008| ) and two Chandra 
observations (|Dai et al.l |2003| ) In this paper we study ten 
archival Chandra observations ranging from January 2006 
until January 2007 and compile the first X-ray light curve 
of Q2237+0305. 
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2 CHANDRA X-RAY DATA 

The data we use was t aken by the Chandra X-ray Observa- 
tory (jWeisskopf et al ]|2P00) in the period between January 
9, 2006 (Juhan date 2453745.0) and January 14, 2007 (Ju- 
lian date 2454114. 9)(PI: Kochanek). The exposure time was 
~ 8 ks per observation. In Table [1] we list the ten obser- 
vation IDs, the exposure times and observation dates. The 
observation IDs are listed in order to reference the individ- 
ual images. We obtained the data from the Chandra Data 
Archive. Figure [T] shows all ten observations. The brightness 
variabilities of the four images are already visible. The four 
quasar images are labeled in the common way from A to D. 



2.1 Chandra configuration 

All ten images were obtained with the ACIS-S detec- 
tor in VFAINT mode. ACIS-S is one of the two de- 
tectors in_the_ACIS (Adva nced CCD Imaging Spectrom- 
eter) (|Garmire et al.l |2003| ) instrument. In our data set 
Q2237-I-0305 was placed on the back-illuminated chip S3 
(1024 X 1024 pixel). 

The VFAINT mode provides the observer with the event 
position in detector coordinates, the event amplitude, the 
arrival time, and the pixel values in a 5x5 island (compared 
to the 3x3 island in FAINT mode). In our case a subarray 
of only half of the ACIS-S array was used for the benefit of 
reducing the read-out time to 1.7 seconds avoiding pile- up 
effects. 



2.2 Data processing 

To process the data we use the CI AO (|Fruscione et al.ll2006l ) 
software package (version 4.2). All ten images are calibrated 
according to the latest Calibration Database (CALDB, ver- 
sion 4.2.0) offered by the Chandra X-ray Center. We re- 
process the images using the standard level 1 event lists 
produced by the Chandra pipeline processing, including the 
latest gain maps and calibration products to create new level 
2 event files. However, we do not apply the VFAINT back- 
ground correction but only the FAINT correction as the 
photons marked as background events in the VFAINT cor- 
rection clearly clump arround the source location, indicating 
that in the 5x5 island, real events of the closely separated 
four quasar PSFs have erroneously been flagged as back- 
ground. In any case, this affected only about ten photons. 
Simulations (PIMMS v3.9) show that for the brightest image 
(image A) in the brightest exposure (6839), (sec Fi gure [p 
pile- up was only 4% using the spectrum bv lDai et alj(|2003l l. 
We also check for contaminating flares during the exposure 
of the single images and find that all ten images show fiat 
light curves. We check the light curves for both, source pho- 
tons (photons in the energy range between 0.5 keV and 8 
keV) as well as background events (photon energies higher 
than 9.5 keV). Finally, we consider an energy range from 
0.5 keV to 8 keV in this paper. The event flies are converted 
into images with a pixel size of 0.492". 



3 PHOTOMETRY 

3.1 The Chandra PSF 

For the purpose of doing photometry on each of the four 
images we created a suitable point spread function (PSF) 
in order to later use it for P SF-fltting. To account for 
Chandra's special Wolter optics (|Wolterll 19521 ) that is com- 
prised in the High Re s olution Mirror Assem b ly (HRMA) 
dWeisskopf et al.l l2000l : ISchwartz et"all l2000l : iGaetz et al.l 
I2OO0I ) and its special imaging characteristics as well as for 
the blended nature of the four quasar images (maximum 
separation of 1.8") that requires a very precise spatial res- 
olution, we create a special PSF. We use the two tools es- 
pecially made for t his scope: ChaRT (Chandra Ray Tracer) 
iCarter et al.ll2003l ). a web interface to a raytrace code de- 
veloped by the Chandra X-ray Center and MARX (Model 
of AXAF Response to X-rays, AXAF w as Chandra's name 
close before launch) (|Wise et al.l 119971 ). ChaRT allows to 
simulate HRMA PSFs at any off-axis angle and for any en- 
ergy or spectrum. 

3.1.1 Creating the PSF 

We use ChaRT to simulate three monochromatic PSFs for 
the off-axis position of Q2237+0305 each with the maximum 
ray density available in order to get a good signal-to-noise 
ratio. We do not simulate pile-up effects. The photon ener- 
gies of the PSFs are 0.5 keV, 2 keV and 8 keV. We choose 
the energies with respect to the spectral range and their rel- 
ative contribution to the spectrum of Q2237-f0305, which 
will be explained more detailed below. The output of ChaRT 
is processed with MARX which projects the rays onto the 
detector taking into account any detector effects. In addi- 
tion to simulating the detector response, MARX uses the 
ray weights to account for the mirror effects, i.e., different 
efficiencies of different shells at different angles and energies. 

3.1.2 Combining the PSFs 

For the PSF fitting we calculate the sum of the three PSFs 
weighting each energy according to its fraction in the spec- 
tra of Q2237-I-0305 as it is seen by Chandra. Th e weight 
of the single PSFs is estimated using XSPEC jA maud 
199(| for a photon index of F = l-90t!5 |^^ (|Dai et all 
20031) and a ga lactic absorption of Nh = 5.5 x lO'^" cm~^ 
(|Dickev fc Lockman 1990 ). The resulting relative weights of 
the PSFs are 1.000 for the 0.5keV PSF, 0.103 for the 2keV 
PSF and 0.008 for the 8keV PSF, respectively. 

The final PSF is dominated by the 0.5 keV PSF, how- 
ever, it has slightly broader wings. It is subsampled by a fac- 
tor of two (corresponding to a binning of ~ 0.25" per pixel) 
and trimmed to a size of ~ 25" x 25" (102 px x 102 px) 
to facilitate the shifting of the PSF. We use this PSF to 
construct the quadruply PSF configuration of our images. 

3.2 PSF fitting 

As the four images of Q2237-1-0305 are very close to each 
other (within a circle of 1.8" diameter), their PSFs are 
blended and we need a sophisticated program to disentangle 
and isolate the individual fiuxes. Thus, we choose GALFIT 
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Figure 1. All ten Chandra images of Q2237+0305 after recalibration and in the energy range between 0.5 keV and 8 keV as used here. 
Each image has its observation ID attached in the upper left. The brightest pixel of image A in January 2006 has 33 counts while the 
brightest pixel in image C has 6 counts. The greyscale is identical for each observation. Brightness variations in image A are clearly 
visible. 



Table 1. The exposure times and starting dates of the Chandra observations of Q2237+0305 used here. 
The observation IDs are listed for future reference. 



Observation ID 


Exposure Time/s 


Start Date 


6831 


7263.84 


2006-01-09 23:39:56 


6832 


7936.11 


2006-05-01 00:46:10 


6833 


7952.71 


2006-05-27 09:45:50 


6834 


7937.77 


2006-06-25 17:24:13 


6835 


7871.38 


2006-07-21 12:03:49 


6836 


7927.81 


2006-08-17 03:59:06 


6837 


7944.41 


2006-09-16 04:50:17 


6838 


7984.25 


2006-10-09 03:41:03 


6839 


7873.04 


2006-11-29 20:08:05 


6840 


7975.95 


2007-01-14 22:30:10 



2.0.3 ijPeng et al.ll2002l ). GALFIT is a two-dimensional fit- 
ting program originally designed to extract structural com- 
ponents from galaxy images. Nevertheless, as GALFIT is 
able to fit several components simultaneously, i.e., quasar 
images, it is an appropriate choice in particular because it is 
able to fit user-provided PSFs. Compared to counting pho- 
tons in circular apertures we can measure the flux of the 
whole PSF and obtain very precise limits of the order of 
0.05"on the position of the PSF on the quasar images. In 
our case, we fit four user-provided PSFs simultaneously in 
order to get the fluxes of each individual image A, B, C and 
D. The four PSFs have fixed relative positions while the ab- 
solute position of the template of the four PSFs is fitted. 
Here, the absolute position is the origin of the relative PSF 
positions listed in Table [21 i.e. the position o f quasar image 
A. We use the relative positions obtained bv lBlanton et all 
l| 19981 ) from U V data because these are minimally affected 
by galaxy light. iBlanton et al.l (|l998l ) determined the astro- 
metric properties of the lens system with the HST WFPC2 
camera in the UV (F336W and F300W filters). We also 



Table 2. Astrometry of the quadruply imaged quasar 
Q2237-I-0305: The table lists th e relative positions of the 
four quasar images as presented by iBlanton et al.l l ll99Sh 



Component 


Right Ascension/" 


Declination/" 


A 


0.000 ±0.0015 


0.000 ±0.0015 


B 


-0.671 ± 0.0015 


1.697 ±0.0015 


C 


0.634 ±0.0015 


1.210 ±0.0015 


D 


-0.866 ± 0.0015 


0.528 ±0.0015 



adapt GALFIT for the X-ray regime to fit fluxes instead 
of magnitudes. 

3.2.1 The fitting algorithm 

The fitting algorithm used in GAL FIT is bas e d on 
the Levenberg-Marquardt algorithm (|Press et al.l Il992l : 
iBevington fc Robinsonir2003l ) that provides a numerical so- 
lution to the problem of minimising a function over a set 
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of nonlinear parameters of the function. By default, GAL- 
FIT and its implemented Levenberg-Marquardt algorithm is 
using a least-squares minimization that is based on the as- 
sumption of normally distributed photon counts. However, 
to cope with the purely Poissonian nature of the X-ray sig- 
nal we fit the data by using a statistic that is based on a 
Poisson distribution. F or this we a dapt GALFIT to mini- 
mize Cash's C-statistic (ICashl l 19791 ) (see Appendix) i nstead 
of the original x^-statistic. Indeed, iHumphrev et al] (|2009| ) 
demonstrate that x^-statistics can give rise to intrinsically 
biased parameter estimates and that the use of Cash's C- 
statistic gives comparatively unbiased parameter estimates. 
In our case we cho ose a sl ightly modified form as imple- 
mented in XSPEC l|Arnaudli 1996.1 1: 



C = -2 ■ ■ log{y^{z,ai, ..,aM)) - yi{z,ai, 

+ yi-yi- log(y») 



, Am ) 



(1) 



Here yi is the i-th datapoint that is fitted to the model 
yi{z, ai, .., au)) with M adjustable parameters aj. The C- 
value is an indicator for the goodness-of-fit and is smallest 
for the best-fitting model. 



3.2.2 The fitting method 

To fit the four PSFs and thereby estimate the fluxes of the 
four images, it is sufficient to choose relatively coarse es- 
timates for the initital fluxes and the position. GALFIT 
then typically converges after a few iterations. It deter- 
mines the PSF template position with an accuracy of of 
a pixel (0.05"). The convergence is robust against changes 
in the initial conditions. Nevertheless, we confirm the fitting 
results by manually stepping through the six-dimensional 
(four fluxes and two coordinates) parameter space. 
Once the best fit is found, we use the fluxes provided by 
GALFIT. The method was carried out on all ten Chandra 
images. Figure [2] shows the raw Chandra image (ID 6838) 
and its best-fitting model. It is notable how well the Ein- 
stein Cross is reproduced and how good the model reflects 
the relative brightness of the four quasar images compared 
to the raw image. 

To estimate the error on the fluxes we manually deter- 
mine the 68% contours for each quasar image flux in each 
epoch using the AC = 1 contour with respect to the mini- 
mum C- value found for the best-flt. 



number. The continuous normal distribution does not have 
this limitation, furthermore the normal distribution even al- 
lows negative values, an aspect that clearly rules out 
based fltting of low number counts. The asymmetric Poisson 
distribution assures that the fitting takes better care of pix- 
els with a low number of counts (< 3 counts) and therefore 
provides an accurate description of the pixels. This is very 
important because all of the quasar image PSFs have pixel 
counts in exactly that count regime. 



4 RESULTS 

The final fluxes for all ten observations and all four lensed 
quasar images in the 0.5 keV - 8 keV band are shown in 
Table [31 The visualization of the data is shown in Figure [3] 
where we present the X-ray light curve that shows the flux 
variations for the ten epochs spread over one year of quasar 
image A (circle), B (diamond), C (square) and D (cross). 

Although the light curve is only sparsely sampled it 
shows flux variabilities in a wide range and for each of the 
four quasar images. Most prominent are the variations in im- 
age A and the general and very steep increase in the fluxes 
of all four images for the last two epochs. As we are look- 
ing at images originating from a single quasar, variations in 
all four images are due to processes in the quasar itself be- 
cause the time-delay is expected to be of the order of a few 
hours only (see below). However, stars in the lensing galaxy 
may also induce brightness variations through microlensing 
which hence are different for each of the quasar images. In 
the following we will separate both types of variability. 



4.1 Microlensing 

Microlensing in the optical is a well known p henomenon 
i n Q2 2 37-I-0305 (e.g. recent st udies including Kochanel3 
(|2004): lAnguita et al.l lISOOSl) : lEigenbrod et al.l (|2008bl h 
iPoindexter fc Kochanekl l|2010l )). Thus we here analyse the 
X-ray light curve to flnd evidence for X-ray microlensing in 
'The Einstein Cross'. 

If part of the flux variations are due to microlensing, the 
pairwise flux ratios of the quasar images should not be con- 
stant. We therefore have to consider that the flux F of each 
image X at time t is a result of different factors. 



X 

macro 



X 
icro 



(t) ■ Q{t) 



(2) 



3.3 Robustness of the fit 

With the implementation of Cash's C-statistic in GALFIT 
we found an ideal way of fltting the X-ray data. The na- 
ture of our signal is clearly Poissonian and this is reflected 
in Cash's statistic. For example, the total number of counts 
for observation 6831 is less than 300 with 33 counts in the 
brightest pixel. 

One of the crucial facts is that the Poisson distribution is 
discrete since the counts in one pixel can only be a natural 

^ Please refer to the XSPEC manual for a more detailed deriva- 
tion. 



Hence the brightness is affected by: 

• Q{t): The time-dependent flux of the quasar as the 
source of the X-ray radiation. It is the same for each quasar 
image. 

• Mmacro- The maguiflcatiou factor that originates from 
the macro lens model and its local lensing parameters k and 
7. It therefore depends on which quasar image we are looking 
at. 

• Mmicro(^)- The Variable maguiflcatiou factor due to mi- 
crolensing changes with time and is different for each quasar 
image. 

We do not need to correct for time-delays be- 
tween the images as they are negligible (|Schneider et al.l 
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Figure 2. Comparison between the central 8x8 pixels in the raw Chandra image (6838) of Q2237+0305 
in the 0.5 keV - 8keV (left panel) and its best fitting model (right panel). The greyscale for both images is 
identical. The crosses mark the positions of the quasar images as obtained in the best fit. 



Table 3. X-ray fluxes of the four images of Q2237-I-0305 derived from the best-fit GALFIT runs. The 
flux for each observation and each lensed quasar image A, B, C and D in the 0.5 keV to 8 keV band is shown. 



Observation ID 


Julian Date 


A 


B 


C 


D 








count rate [10 * • s ^] 




6831 


2453745.0 


170.7 ± 17.9 


129.4 ± 16.5 


38.6 ±8.3 


57.8 ± 13.8 


6832 


2453856.0 


331.4 ±23.9 


157.5 ± 16.4 


120.1 ± 15.1 


92.0 ± 15.1 


6833 


2453882.4 


182.3 ± 16.4 


78.0 ± 11.3 


66.6 ± 11.3 


36.5 ± 10.1 


6834 


2453911.7 


405.7 ±25.2 


162.5 ± 17.6 


118.4 ± 15.1 


90.7 ± 16.4 


6835 


2453937.5 


471.3 ±26.7 


115.6 ± 14.0 


80.0 ± 12.7 


85.1 ± 15.3 


6836 


2453964.2 


273.7 ±20.2 


85.8 ± 12.6 


70.6 ±11.35 


59.3 ± 13.9 


6837 


2453994.2 


244.2 ± 22.6 


102.0 ± 12.6 


57.9 ± 11.3 


61.7 ± 12.6 


6838 


2454017.2 


231.7 ± 18.8 


85.2 ± 12.5 


72.6 ± 11.3 


43.8 ± 12.5 


6839 


2454068.8 


882.8 ±36.8 


323.9 ±22.9 


176.5 ± 19.1 


165.1 ± 22.9 


6840 


2454114.9 


670.7 ±36.4 


214.4 ±20.1 


185.6 ± 18.8 


131.6 ± 21.3 



19881: iRixetaLl Il992l; IWambseanss fc Paczvnskil 1 19941 : crolensing signal in image D when compared with the ratio 



Schmidt et aIlll998l : [Dai et al.ll2003l : iKoptelova et al.ll2006l ) 



As we are interested in microlensing we form the pairwise 
flux ratios between all quasar images X and Y {X, Y — 
1,..,4; X^Y): 



^ (j'^ /^macro ' /^micro(^) ' Q(^) /^macro ' /^micro(0 

Mmacro • Mmicro(*) ■ Mmacro ' Mmicro (*) 



(3) 



This way the intrinsic fluctuations of the quasar cancel 
out. The results are shown in Figure |4] where we logarith- 
mically plotted the pairwise flux ratios of the quasar images 
for every observation and according to equation ([3]) . The dif- 
ference to a light curve including the macro magnification 
Mmacro IS just a coustaut offsct iu the logarithmic plot and 
is of no importance for our analysis. Additionally, we also 
plotted the MIR flux ratios (Agol et al. 2000), i.e. the macro 
magnification ratios, for each pair of quasar images. 

In Figure|4]we see that the ratios and ^ show a sim- 
ilar variation signature (with the exception of the first data 
point), i.e. an increasing flux ratio that has a maximum in 
the flth epoch and decreases afterwards. This behaviour is 
not mirrored by the ratio ^ for these epochs, which supports 
the view that the maximum in the ratios and ^ is caused 
by microlensing in image A. The variations of ratios and 
S are consistent with an independent low-amplitude mi- 



^. However, the evidence for microlensing is most promi- 
nent in image A and the ratio which is why we only focus 
on this quasar image and image ratio, respectively. 

The arguments given above and the fact that the cor- 
related variations especially in the ratio do not support 
random fluctuations as the cause, leads to the conclusion 
that we have observed microlensing induced variabilities. 
However, by comparing the observed flux ratios with the 
macro magnification ratios, we find that all quasar images 
may be affected by microlensirig because their ratios are not 
consistent with the MIR dataQ The fiux ratios can also be 
affected by substructures in the lensing halo, however, due 
to the higher mass scale the time-scale of the variations in- 
duced by substructures is much longer (due to the higher 
mass-scal e) than the time-sca le considered in our analysis 
(Mctcafl 2 0051 : IXu et"al]|2010t ). 

Focussing on th e ratio of image A to im a ge B and com- 
parin g it to OGLE (|Wozniak et all l2000bl lal: lUdalski et al.l 
I2OO6I ) observations (V-band, restframe UV) obtained at the 
same time, we find that the maximum as well as the subse- 
quent flattening of the ratio in the X-ray regime coincides 



^ Note that the offset between the X-r ay and IR fluxes can b e 
slightly offset by differential absorption jEigenbrod et al.ll2008al) . 
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Figure 3. The final X-ray light curve of the four images of Q2237+0305 in the 0.5 keV - 8 keV band as obtained with the C- 
statistic based GALFIT. The figure shows the variations with time for all four lensed images. The datapoints and errorbars have 
been connected with shaded bands for each quasar image. A sudden increase of flux due to intrinsic variations in 'the central engine' 
of the quasar can be seen for the last two epochs. The corresponding data can be found in Table |3] 



with the optical observations. This is shown in Figure [5] 
where the ratio in the X-rays is overplotted with the 
OGLE data. For the most prominent peak around Julian 
date 2453937.5, the increase since Julian date 2453856.0 of 
the X-ray ratio A/B is larger by a factor ~ 1.6 compared to 
the signal in the optical. Furthermore, also the A-B differ- 
ence light curve in the UV (Eigcnbrod ct al. 2008b) for that 
period shows similar characteristics. However, the events are 
less distinct at UV and optical wavelengths which is exactly 
what we expect as the UV and the optical emitting regions 
are thought to be larg er than t he X-ray emi tting region 
llWambs g anss fc Paczvnski 1991i: iRauch fc Blan dford 199L; 
iDai et al.l [20031 : iPoolev et al.1 12OO6I : iKochanek et al.1 \26m . 
The inset in Figure [S] shows theoretical light curves for two 
Gaussian sources crossing a straight caustic. The radii are 
different by a factor of 6.7. 



4.2 Intrinsic variations 

The fact that there is strong evidence for microlensing in im- 
age A and D does not explain the high flux rates for the last 
two epochs and the flux rate drop for the third observation 
as the microlensing in image A had its peak at JD 2453937.5 
(July 2006) . These fluctuations can only be explained by in- 



trinsic variations of the quasar itself. X-ray variability on 
such short time scales of weeks indicates that the emission 
generation takes place in a relatively compact region and 
that the emission processes, that make quasars so luminous, 
are highly variable. Such strong variations are not seen in 
the optical OGLE light curve. Merely, the minimum at JD 
~ 2454010 can be seen. 



5 SUMMARY AND OUTLOOK 

In this paper we have analysed archival Chandra data of the 
gravitationally lensed quasar Q2237-I-0305. We compiled an 
X-ray light curve for all four images. The data set comprises 
ten epochs ranging from January 2006 to January 2007. 
Because of the blended nature of the four images (see Figure 
[1]) it was necessary to simultaneously fit appropriate PSFs 
to all four images in order to obtain proper photometry. 
For this, we used a simulated PSF which accounts for 
the optical properties of the Chandra observatory and the 
spectrum of the quasar X-ray emission. The fitting was 
accomplished by a two-dimensional fitting algorithm and 
a grid search by minimising Cash's C-statistic. Finally we 
analysed the light curve and found evidence for microlensing 
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Figure 4. Pairwise X-ray count rate ratios of the four quasar images A through D of Q2237+0305 in order 
to identify microle nsing. The rati os are calculated according to Equation [3] The liorizontal lines mark the 
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variations in quasar image A (see Figure 2]). The X-ray 
microlensing signal in image A coincides with the signal in 
the optical OGLE light curve. Assuming that this parallel 
behaviour is caused by the same process, i.e. one source 
that is microlensed, the amplitude of the microlensing 
signal is a direct measure for the source size in the 
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1991 




2008 



Eigenbrod et al.1 l2008bl )~ As the microlensing signal in the 



X-ray regime is much more prominent than in the optical, 
this suggests that the X-ray emission region is much smaller 
than t he optical emission region of the quasar (jPoolev et al.l 
I2OO9I : iDai et alJlioiol ). In a future paper (Zimmer et al. 
2010, in prep.) we will make use of this effect to measure 
the size of the X-ray and the optical emission region in 
Q2237-I-0305. While the time-delays in Q2237-I-0305 are 
negligible the intrinsic variations have to be considered 
in the analysis. Figure [3] clearly shows how strong these 
variations are in the X-ray regime. We do not find such 
strong fluctuations in the optical. This indicates that the 
mechanism powering the quasar is variable and leads to 
brightness variations by a factor of ~ 4 on time-scales of 
less than 30 days. 
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p{y,) = i/-l^,ffli.-,"Mj e ^^^^ 



y^'■ 



So the probability of the whole data set is the product 
over all datapoint probabilities p{yi). 



(A2) 



To maximize the data set probability by estimating 
the model parameters oi, ..,aM one can also minimize —2 • 
log(P): 



-2 ■ log{P) = -2 ■ YJJ^ ■ log(yi(z, ffli, ••, flAf))- ^^g^ 

yt{z,ai, ..,aM) - log(2/i!) 
leading to the definition of C: 



C = -2-^y» •log(j/i(z,ai, .., au)) ~ yi{z, ai, ..,aM) (A4) 

i = l 

The term — log(yi!) can be left out as it does not depend 
on the model parameters and is therefore just a constant off- 
set. For our implementation in GALFIT we chose a slig htly 
modified form as implemented in Xspec (|Arnaudlll9m 



C = -2 ■ ^yi-\og{yi{z,ai, ..,aA/)) - yi{^ 

+ yi-yi- log(j/i) 



z, ai, .., aM 1 



(A5) 



This translates to the following function to be mini- 
mized by GALFIT: 



C = —2 ■ ^ y~^fiux3:,y • log(modeL,j;) — modeL,;/-!- 

1 = 1 y = l 

fiUXj;,y - HUXx^y ' lOg ( fiUX^; , y ) 



This paper has been typeset from a T^j]X/ ffl^jX file prepared 
by the author. 



APPENDIX A: CASH'S C-STATISTIC 

In comparison to x^-statistic Cash's C-statistic l)Cashlll979l l 
is bcised on the assumption that each datapoint j/i is a ran- 
dom variable out of a Poissonian distribution with an ex- 
pected value of yi{z, ai, .., om)- Therefore, it has the proba- 
bility: 
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